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SUMMARY 


The report presents information on the stress problems in the 
analysis of pressurized cabins of high-altitude aircraft not met with 
in other fields of stress analysis relating to aircraft. The material 
may be roughly divided into shell problems and plate problems, the 
former being concerned with the curved walls of the cabin or pressure 
vessel and the latter being concerned with small rectangular panels of 
its walls, framed by stiffeners, but not necessarily plane. 


INTRODUCTION 


The analysis of pressurized cabins of high-altitude aircraft pre- 
sents particular stress problems, not usually met with in other fields 
of stress analysis relating to aircraft. It is the purpose of the 
present report to gather information on these problems and to make it 
easily accessible to aircraft engineers. Some of the work in this field 
is presented in references 1 to 10. 

This report contains a choice of subjects taken from the theory of 
plates and shells which is expected to be useful for the designer of 
pressurized airplane cabins or similar lightweight pressure vessels. 

This material may be roughly divided into shell problems and plate ' 
problems, the former being concerned with the curved walls of the cabin 
or pressure vessel and the latter, with small rectangular panels of its 
walls, framed by stiffeners, but not necessarily plane. 

As far as shell problems are concerned, some use has been made of 
a manuscript for a book on "Stresses in Shells," which the author is 
preparing. (See reference 3.) The prospect that this book will be 
available some time in 1952 makes it possible to discuss in the present 
report several problems which are too complex to explain here in all 
their mathematical details. 

The pressurized cabin is a rather new element in the airplane 
structure and will, in all probability,- undergo future development. In 
view of this situation, no attempt has been made to present anything 
like a textbook on the subject giving time-tested methods for solving 
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various problems, but rather an attempt has been made to show the general 
lines of thought which have proved to be useful and to give suggestions 
for their application. 

This investigation was carried out at Stanford University under the 
sponsorship and with the financial assistance of the National Advisory 
Committee for Aeronautics. 


SYMBOLS 


x,y>z 

0,0 

u,v,w 

a 

a,b 

l 

t 

P 


X,Y,Z 


W N x 





or 

T 


rectangular coordinates 
angular coordinates 

displacements ' t 

radius of cylinder or sphere 

sides of rectangle; axes of ellipse or ellipsoid 
span of beam 

thickness of plate or shell 

pressure difference between interior and exterior of cabin 

distributed load on shells (force per unit area of middle 
surface), in directions 0, 0, and radial 

normal forces in shells (force per -unit length of section), 
in direction 0, 0, or x 

bending moment in plates and shells (moment per unit 
length of section) 

twisting moment in plates (moment per unit length of 
section) 

normal stress 

shear stress 


E elastic modulus 

V Poisson' s ratio 
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SHELL PROBLEMS 
Cylindrical Shell 


Circular cylinder .- The fuselage of a high-altitude passenger plane 
is usually of circular cross section and is, for most of its length, 
almost cylindrical. Some useful information regarding its strength may 
he found, therefore, when a circular cylinder closed at both ends by 
some kind of bulkhead which permits the air pressure inside te be greater 
than that outside (fig. l) is considered. The pressure difference will 
be called p. 

For a homogeneous shell of thickness t the stresses produced by 
this pressure are given by the well-known boiler formulas for hoop stress 
CM and axial stress cr x : 


00 = pa/t 
a x = pa/2t 


( 1 ) 


The shell of a pressure cabin is reinforced by rings and stringers, 
which participate in carrying the load. The stringers will always be 
spaced closely enough to make the distribution of the longitudinal stress 
on the skin between them practically uniform. With the rings this may 
be different. The limiting case, that is, that they too are closely 
, spaced, will be considered here. 

In finding the stresses, start from the internal forces per unit 
length of section acting in the shell. When a slice of length Ax = 1 
is cut out of the shell (fig. 2), the hoop force 

N 0 “ P a 

is found, and when the force pna^ acting on each bulkhead is distributed 
over the circumference 2jta of the cylinder, the longitudinal force 



transmitted by the unit length of a section right across the shell is 
found. 
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If the shell has no stiffeners, the stresses a j and a x are 
found by dividing N0 and N x by the vail thickness t, which, of 

course, results in the boiler formulas (l). In the cabin shell are 
rings of cross section A E at distance 2 from each other and stringers 
of cross section A^ at an angular distance 6 (see fig. 3 ). If these 
areas are distributed over the cross section of the skin, the effective 
thicknesses 


H 


= t 



( 2 ) 


are introduced; however, the stresses 00 and ct x are not simply the 
quotients N 0 /t 0 and N x /t x (see, e.g., reference l). The reason for 

this is apparent when one considers the fact that the skin is in a two- 
dimensional state of stress and therefore for the same strain its stress 
is different from that in the stiffeners. 

Let the stresses in the skin be 00 and a x as before, in the 
stringers, a^, and in the rings, a R . Then Hooke's law will yield the 
following relations for the hoop strain «0 and the longitudinal 
strain * x : 


Ee 0 = CT0 _ V0 X 

= a R 


(3a) 


Ee x = a x - Va0 
" a L 

E being Young's modulus and V Poisson's ratio. 

On the other hand, the definition of the internal forces is: 

Ar 

N0 = tcx0 + -f 0 R 

N x “ He + ^ 


(3b) 
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Solving the four equations (3) and (^) for the stresses, 

+ v ( t 0 - 

^ (l - V 2 )t 0 t x + V 2 t(t 0 + t x - t) 

pa 2 t x + V (t0 - t) 

2 (l “ V 2 )t0t x + V 2 t(t0 + t x - t) 

[u - v 2 )t x + v 2 t] “ v ™ x 

= — 

(l - V 2 )t0t x + + t x - t) 

pa 2(1 - v 2 )t x - v(l - 2v)t 

2 (l - V 2 )t0t x + V 2 t(t0 + t x - t) 


( 5 a) 


( 5 b) 



* v ( t x - 

(l - V 2 )t0t x + V 2 t(t0 


pa 

~ 


t0 + 2v (t x - t) 

(l - V 2 )t0t x + V 2 t(t0 + t x - t) 


( 5 q) 


[(l - V 2 )t0 + V 2 t]w x - VtN 0 
(l - V 2 )t0t x + V 2 t(t0 + t x - t) 


pa I 

(1 

- v 2 ; 

^ t0 - v(2 - v)t 

2 (1 -V 2 ) 


0 

+ V^t | 

(t0 + t x - t) 


( 5 d) 
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When the rings are far apart, these formulas are no longer appli- 
cable * The problem must then be split, with the shell without rings 
considered first and the influence of the rings introduced afterward 
(see section entitled "Interaction between Shell and Rings" ) . When 
there are no rings t0 = t, and the formulas are simplified considerably: 



(6a) 


pa(l - 2v) Vpa 



(6b) 


°L 


(l - 2v)pa 

2t x 


(6c) 


It appears that op, is always smaller by a factor 1 - 2v than it 
would be if it were obtained by simply distributing N x over the whole 
section. For the skin stress cr x the factor depends on the ratio 
A]Ja5t, and if one writes 



the factor k will be as shown in figure 4. For AL/a&t = 0 the 
boiler formulas are valid, and a x = O.5a0. For Ap^abt = 1.0, the 
diagram shows cr x = O.4cr0. The difference between these two values 
of cr x is small, but both are much less than the hoop stress. This 

is very desirable since the over-all bending of the fuselage due to air 
forces acting on the control surfaces produces additional stresses a x 

which must be superimposed on the stress cr x due to cabin pressure. • 


Since the stringers take an important share of the axial load, it 
is not good practice to interrupt them at the rings. Care should be 
taken to insure that the forces carried by the stringers can go straight 
through from one bulkhead to the other, or to the end of the cabin shell. 

Double cylinders .- The circular cross section is certainly the best 
one, both for aerodynamic and structural reasons. However, it has some 
practical disadvantages when used as a. passenger cabin. Most serious 



NACA TN 2612 


7 


is the fact that a horizontal floor must necessarily he built in, 
requiring additional weight and leaving beneath it space which is not 
easily used. 

This situation is improved by a cross section which, with some 
exaggeration, is shown in figure 5. It consists of parts of two circles 
and a straight horizontal tie between them. 

Begin with a discussion of the weight of this structure. Under the 
action of an internal pressure p the hoop stresses in the upper cylin- 
der ow.. and in the lower cylinder ctj 0 will be: 


^1 - 

*02 " 2 * 2 1 *2 

The stress in the tie follows from the equilibrium at its ends (fig. 6): 


ff 3^3 = a 01^1 cos a l + ff 02^2 cos ^2 
= cos a l + a 2 cos a 2 ) 


If t-p and t^ are chosen such that the three stresses are all 

equal to the same value a given by the allowable stress in the material. 




cos an + &o cos 


a 2 ) 


The material invested in the structure is given by the area A^^ 
of its metallic cross section. The two circles contribute to it 
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2ti a i(rt - a-jJ + 2t 2 a 2 (jt - a 2 ) = glia-^ 2 ^ - a^) + 2a 2 2 (jt - 


and the tie contributes 


2toC 


= 2 

a\ 


a 2 )' 


Now 


aj cos + a 2 cos a 0 )c 


c = a-j_ sin 


= a 2 sin 0-2 


and hence 


2t^c = 2 cos a^_ sin a-j. + a 2 2 cos 02 sin 


■U 


a i sin 2a n + an sin 


202) 


Summing up the three parts , the total metallic cross section is found 
to be: 


A m = 2 


^ja-]_ 2 ^n - a^_ + i sin 2a^ + a 2 2 ^jt - 02 + \ sin ^■ a 2^ 


On the other hand, the area of the hollow cross section A^ 
describes the useful space in the fuselage. It is 

Ah = a-^At - a-j_ + i sin 2a-^ + a 2 ^jt - 02 + ^ sin 2a 2 j 


It is seen that the ratio of the two areas 


^m/A h -2g 
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does not depend upon the particular choice of the dimensions of the 
cross section. For a simple circular cylinder of radius a there is 
obtained by similar reasoning: 



Am = 2 * at 

= 2« — a 2 
(T 

A h = « a2 

and hence 

Am/A h = 2 E 

as before. This indicates that for the same inside space the same 
structural weight is required and one is freed from weight considera- 
tions when choosing that combination of a^, a a^, and a ^ which 

seems best for other reasons The validity of this result is restricted 
to cross sections where the tie acts in tension, and this is exactly the 
configuration which is most interesting in aircraft construction. In 
practical applications, of course, additional stresses will change the 
picture to some extent, and the weight of different shapes will not be 
equal, but the important fact remains that there is no first-order loss 
or gain in choosing one or another of the sections compared. 


Interaction between Shell and Rings 

Bending of a cylindrical shell .- If the rings are not spaced closely 
enough to be considered as part of an anisotropic shell, the problem 
illustrated by figure 7 must be treated. Cut the shell in the plane of 
the ring and at its connection with the ring. The pressure p applied 
to the shell will lead to a hoop strain which may be found from 

Hooke's law (3) and formulas (6a) to (6b) and consequently will lead to 
a radial displacement 


..1 


/ 


w o “ ae 0 

1 - 2V) \ 
2t x ) 


pa 2 (l - V 2 v( 
~\ t ’ " ~ 
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The ring receives no load and, therefore, has no deformation. In order 
to close the gap between the ring and the deformed shell, add shearing 
forces T per unit length of the edges. In the ring of cross section 
Ar these shearing forces produce the stress 


2Ta 



and hence the radial displacement 


W R 


2Ta 2 

EAr 


( 7 ) 


p 

(More exactly, aa-^ should be written instead of a c , where is the 

radius of the center line of the ring. ) 


For the shell, the force T is a transverse shear Q x which pro- 
duces bending stresses. In order to find them, some details of the 
theory of bending of an anisotropic cylinder must be developed. It is 
necessary to consider only the internal forces and moments shown on the 
shell element in figure 8: The hoop force N^,- the bending moment M x , 

and the transverse shear Q x . They are all functions of x (fig. 7)* 
as is the radial deflection w. 


The forces and moments must satisfy the conditions of equilibrium 
of the shell element. They yield two equations: 



(8a) 


d ^x 

dx + ^0 = 0 

which, after elimination of give the relation: 


(8b) 


d 2 M,, . 

a + Ngj = 0 (9) 

dx 2 

The hoop force produces a hoop strain S0 which may be obtained from 

equations (3a) and (5a) to (5c) with N x = 0. This strain leads to a 
radial displacement 
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w = ae 


0 


= N 0 a 
H0a 
= D 0 


t x - v 2 (t x - t) 


Ett^ 


t x 

The constant B0 = Et 5— -, which has the dimension of a force 

t x “ v (^x " ’ t ) 

per -unit length, is the extensional rigidity of the shell in the direc- 
tion of the hoop forces. Figure 9 shows that, in the range of practical 
interest, D0 is only slightly greater than Et, and it is safe to say 


N0 = Et — , (10) 

The bending moment M x produces a curvature d 2 w/dx 2 of the generators. 
If I is the moment of inertia of the cross section of a stringer and 
the attached skin of width aS divided by the distance a5 of the 
stringers, . . 

= El ~ (11) 

dx 2 



Here too the coefficient El is slightly influenced by the fact that the 
skin has a two-dimensional stress system. This refinement of the theory 
will not be discussed here. There is another circumstance, perhaps even 
more serious, which will also be neglected here: The centroid of the 

section to which I is referred is not exactly at the distance a from 
the axis of the cylinder but at a somewhat shorter distance. This 
influence may be studied with a more general set of equations, but since 
the difference of the two radii is not great, it will probably not be 
of first-order importance: however, it may be responsible for some 

second-order effects which otherwise might not be explained. 

By introducing the expressions for N0 and into equation (9), 

the differential equation of the problem is obtained: 


- El 


d**v 


dx 


+ Et 


w 


0 


( 12 ) 
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The general solution of this equation consists of four terms. Only 
those which are symmetrical with respect to the plane x = 0 are 
needed. They are: 


KX KX 


kx . KX 


w = Ct cosh — cos — + Co sinh — sin 

-L a a -c a a 


with 


k 

■-N 


ta 

41 


The boundary conditions at x = Z/2 are that the slope dw/dx must be 
zero an d that the deflection must assume a certain value w^. (This 

will be discussed later. ) 

(**) = 0 

Wx-l/2 


v x=2/2 = W 1 


Introducing the solution here, Cq and C2 are found and then 


wi 


w = 


cosh p sinh p + cos p sin p 


0 


cosh P sin p + 


sinh p cos p) cosh cos + ( cosl1 P sin P 


. , n 0 n . , 2px . 2px 

sinh p cos P) smh sxn 

l l 


P = 


kZ_ 

2a 


I >/t 


2^ VI 


with 
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The internal forces can now be found easily. From equations (8a) and (ll) 
it follows that 


Qx = 


^x 

dx 


= El 


d^w 

dx3 


EtZw-i 


2a%(cosh p sinh 3 + cos p sin p) 


( . . . _ . 23 x . 2px 

cosh P sin P cosh — — sin — — + 


2Bx 

sinh p cos p sinh cos 


¥) 


For x = l/2 this is the force T applied to the edge of the shell: 


T 


EtZwi 


cosh 2 p - cos 2 P 


2a 2 p cosh p sinh p + cos p sin P 


(13) 


Combining this with the preceding formula 


Ox 


= T 


cosh p sin p cosh s ^- n + P cos 0 sinh cos 

cosh 2 p - cos 2 p 


Introducing the solution w into equations (ll) and (10) and expressing 
w 1 by T, • 


M x = El 


d^w 

dx^ 


= -T 7— - 


- — — R 

- cos^-B L. 


cosh p sin p 


^P cosh 2 p - cos 2 P 
inh P cos p) cosh cos ^2. 


(cosh p sin p + sinh p cos p) sinh sin 

z *• 
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and, adding the hoop force due to the pressure p, 


Nwf = pa + 


Etw 


= pa - T 


2ag_ 1 

1 cosh 2 3 - cos 2 3 


J^cosh 


p sin 3 + 


. . „ , 23x 23x , 

sinh 3 cos 3) cosh -j- cos — + 


(cosh 3 sin 3 - sinh 3 cos 3) sinh 


£3x 

l 


sin 



The magnitude of the shearing force T still has to he found. Write, 
for abbreviation, T = -kw-j_, where k is defined by equation (13). 

Then the following deformations are found: Under the action of the 

internal pressure alone the shell has the deflection w 0 given at the 
beginning of this section, and the ring, none. The additional load T 
bends the edge of the shell back, producing w-j_ = -T/k, and the ring 

has a positive displacement w R = 2Ta^/EA R as was seen earlier. 


Now, under the combined action of pressure p and shear T, the 
radial displacement of shell and ring must be -the same 



2Ta^ 


2a 2 + l 
EA r k 


From this equation 
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Upon the introduction of w Q 


and k this yields 



p /l - v 2 y(l - 2 v) \ 

2\ t 2t x ) 

JL ( c °sh p sinh p + cos p sin p 
\ cosh 2 p - cos 2 p 


(l4) 


With the numerical value of T from this formula, one may obtain from 
the preceding formulas values for and M J{ . The complete solution 

of the mechanical problem of the interaction of the shell and the rings 
is now obtained. 

The shear depends on many parameters, and no attempt has been made 
here to represent formula (l4) by diagrams. However, the distribution 
of N 0 and ^ along a generator of the cylinder depends essentially 

on p. When P is small (closely spaced rings, heavy stringers), the 
picture looks somewhat like figure 10 (a), and the case where it is 
adequate to represent the influence of the stringers by the effective 
thickness t f as defined by equation (2) is approached. But when p 
is great (rings far apart, light stringers), the internal forces are 
like those sketched in figure 10(b): In this case the influence of each 

ring is locally restrained. 

Floating skin .- From figure 10 it appears that there is not much 
virtue in providing rings to help the skin carry the cabin pressure 
because the skin alone can do that well enough and the rings only cause 
trouble. The rings disturb the simple stress system considerably, and 
the force 2T transmitted from the shell to the ring produces a highly 
undesirable tensile stress in the rivets which connect the skin to the 
ring. - ' 

However, the rings are needed for many important purposes. They 
help to introduce the local load gently into the shell, they support 
the stringers against buckling, and they stiffen the shell as a whole 
to prevent a collapse by large-scale buckling. The problem is, how 
does one make the rings available for all these purposes without 
introducing the forces T? „ 

The solution is the floating skin. Its basic idea may be explained 
from equation (l4). If the rings are very weak, Ar 0, the denomi- 

nator of this formula becomes infinite, and T = 0. The term 2/Ar 
comes from vr, equation ( 7 )> and represents the deformability of the 
ring by a radial load T, or, more exactly, the radial displacement of 
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those points where the ring is fixed to the shell (skin and stringers). 
This deformability may be increased easily without weakening the ring. 
It is only necessary to intersperse an elastic element between the ring 
proper and the shell (including stringers) as it is indicated by the 
sketch, figure 11. Formula (7) must then be replaced by 


w R 


2Ta 2 

EA R 


+ w R * 


= ==^- + cT 
EAr 

\ 

c being an elastic constant depending on the shape and size of the 
connecting link between ring and skin. In formula (l4) it is then 
necessary to write 

_2_ Ec 

' a R a 2 


instead of 2/Ar and now there is the possibility of making the denomi- 
nator as great as desired. Of course, such a flexible connection is 
only worth while if the denominator in equation (l^) is appreciably 
increased by adding the term Ec/a 2 . 


Doors and Windows 

When the cabin is under pressure, the door must, of course, be 
closed, but it cannot be expected that the door panel will be very 
efficient in transmitting hoop forces N0 or longitudinal forces N x 

across the door opening. Both have to be carried around it by the door 
frame, and this disturbance of the smooth flow of forces will certainly 
lead to an increase in structural weight. In order to keep this increase 
as small as possible, some details of the local stress system will have 
to be studied. 

Since the door needs a frame, it is reasonable to extend the lateral 
parts of this frame all around the shell as two of its rings. Outside of 
the part of the fuselage limited by these rings the hoop force does not 
meet with any obstruction. The problem is, what must be done with the 
forces which are intercepted by the sill and the head of the frame? With 
the usual dimensions of fuselages and doors these forces are considerable. 
A door frame strong enough to resist them would be a heavy structure. 
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but worse than that, it would not accept the load. It would deflect in 
the direction of the pull, and the deflection would lead to a decrease 
of the pull. 

It seems wise, therefore, to allow the shell itself to do that which 
it can do so easily and to give the horizontal members of the door frame 
only that stiffness which is required to press the door firmly against 
it, that is, bending stiffness against radial forces. It is necessary, 
then, to solve the following mechanical problem (fig. 12): A cylindrical 

shell extending over an angle a < 3^0° is limited by two circular rings 
and by two straight end members. These end members have no rigidity 
against bending in the tangential plane to the cylinder, but they have 
enough cross section to be considered as inextensible for our purposes. 

The shell is subjected to an internal pressure p. 

The stress system set up under these conditions may be split into 
two parts: One is a hoop force N0 = pa, acting everywhere (also on the 

straight edge) and resisting to the load p; the other one is a system of 
internal forces produced by an external load F0 = -pa applied to the 

edge members and canceling there the force pa of the elementary solution. 

The task which is now to be done is to find this second stress system. 
In the theory of shells it is shown that the tangential load -pa cannot 
be carried by the shell without resorting to bending stresses. There are 
different methods of treating this bending problem; this is a simple one. 
Although its application in this case may not be entirely legitimate, 
it will give a fair idea of what is happening, and that should be enough. 

The problem may be reduced to a differential equation for the 
bending moment M0 (for the notations, see figs. 8 and 12; for more 
details see reference 2, p. 139* or reference 3): 


M0 :::: + (2 + v)M 0 " s:: + 2M0 ::: + (l + 2v)?y" ,:: + 2(2 + v)M0 ,,:t + 

,2 


M0 :: + VM0’ M,,, °- + (1 + v)^y m -- + (2 + V)M 0 "‘* + M0 1111 = 0 (15) 


Here dots indicate derivatives with respect to 0, primes indicate 
derivatives with respect to the dimensionless coordinate x/a, and k 
is the shell parameter 


k = 


t 


2 


12a 


2 
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It may easily be verified that 

M0 = Ce^ sin M 

is a solution of the differential equation. Here X is still an 
arbitrary constant. Write 

X = 

l 

where n is a positive integer; the discussion will later be confined 
to n = 1. 

When the solution is introduced into the differential equation, it 
is found that m must satisfy a certain algebraic equation. After some 
drastic simplifications (reference 4) it may be brought into the 
following form: 

' i 

^ . l) 2 ^ + 1 -...y. 2 X 1 * - 0 

k 


This equation has the complex solutions 


with 



i • \ ^(1 - v 2 ) 


(16) 


When any one of the eight complex values m is introduced into the 
formula for M^, one elementary solution is obtained. They all show a 
variability in x-direction according to sin n nx/Z. The same factor will 
appear in the corresponding hoop forces N0. When it is desired to use 
these solutions to describe the stresses due to a uniform distribution 


= -P a 
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of the hoop force on the edge 0 = 0, this distribution has to be 
resolved into its harmonic components: 


-pa = 


%)a 

rt 


itx , 1 . 3 «x 1 . 5«x , 

sm — + — sm + — sin - — + 

13 l 5 Z 




The discussion will be confined here to the first term of this series 
which will show^he essential features of the stress pattern. Corre- 
spondingly, set n = 1. 


The set of eight elementary solutions which is obtained from the 
eight possible, values of m may be replaced by an equivalent set of 
eight linear combinations, each of which is a product of an exponential 
and a trigonometric function of k-j _0 or where 


*1 = £ \J"n cos 22.5° 


*2 = £ \fn sin 22.5° 

Using a suitable set of boundary conditions to determine the constants C 
with which these solutions must be multiplied, solutions for many cases 
of loading and supporting the edge may be found. The full expression 
for the bending moment Mq, and the values of the displacements u (in 

x-direction) and w (radially outward) for the edge cp = 0 are given 
here for three important cases: . 


(l) Normal forces N0 = sin applied to a free (unsupported) 


edge. 


M0 = - — i— ^\|2 + 1^ |e (cos )<20 + sin K 2$) ~ e ^ cos K 10 


• k 2 0 


sm 


itx 

l 

(17a) 


with the displacements at the edge 0=0, 


u = 


EtX 3 


2 + l)\|2£ 2 - 1 - vX 2 cos 22 
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Fja 

W = r; 

EtA. 4 


({2 + l) 2 ^ + (l - 2X 2 )({2 + l))[2C 2 + (2 + V - A 2 )*. 2 ] sin 


(2) Shearing forces N0 X = F2 cos — , applied to a free edge 


F2aX 


*¥= 25 3 


^2 + ^2 je Kl ^ jcos (C20 + (^2 + l) sin k 2 $\ 


e * 2 ^ jcos *^0 + ^\|2” - l) sin *10^ sin ^ 
with the displacements at the edge 0=0, 


Fp a £ 


u = 


EtX' 


5 jzp * W =°* f 


W = - 


F 2 a £ 

EtX 3 


^|2,2 + |2 ^ 2 + 1 - 2X 2 ^ sin 


nx 

Z 


(3) Combined action of a shearing force N0 X = F3 cos — and e 
transverse force = -XF3 sin — , applied to a free edge 


F 3 aX 


fTji jr 15 ^ [(js + l) cos *20 + (2^2* + l) sin. *g0j 


"* ’ 2£ 

e “ K 20 + l) cos *20 - sin K ]_0]^j' 

with the displacements at the edge 0=0, 


JtX 

sm -y 


u = - 


f 3^~ 

EtX 2 


+ \|2 


cos 


nx 


w = 


f 3°^ 3 

EtX 3 


(2^2 + ^2 (J2+ l)? 


' 2 + 1 - 2 X . 2 ‘ -- * X 


2l • nx 
sin — 


nx 

Z 


(17b) 


(17c) 
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Coming back to the original problem represented by figure 12, the 
bending moment M 0 may be found by superposing the three solutions 
(equations (17&) to (17c)) with appropriate values of the constants Fp, 
F£> and F^. These can be found from three boundary conditions at the 
edge 0 = 0 . 

The first harmonic of the load pa shown in figure 12 is 
N0 = - sin — . In order to give N 0 this value, it is necessary 

to set 


Fi = 


4pa. 

Tt 


and this is the first of three conditions. The other two follow from 
the deformation which the door frame imposes on the shell. Since the 
cross section of this frame is assumed to be large enough to neglect 
axial deformations, u = 0 at 0=0 for the shell. Another assump- 
tion formulated previously is that the door frame will not deflect 
very much in the w-direction. Therefore, w = 0 at 0=0 for the 
shell. 

When u and w are expressed as sums of the contributions of Fp, 
1 ? 2 > and. F 3 , according to the formulas given before, and set equal to 
zero, there are two linear equations from which F 2 and F 3 may be 
found. The result is as follows: 


f 2 = 


pa 


^ - (2 (\f 2 + l )£ 4 - (<2 + l)(l + v\ 2 )£‘ 


3 a, 

1 - 2vX 2 + (1 + 2 v)xA] 


F 3 = 


Pa 


^2 - \[2 ft[2 + l) ^ - (l - vA 2 )^ 2 + 1 + 2vX. 2 - (l + 2v)\ k 


The bending moment M0 due to the combined action of Fp, F^, 
and F 3 can now be computed; going back through details of tpe theory 
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which have not been reproduced here, the hoop force N0 and the longi- 
tudinal force N x in the shell can also be found. 

This has been done for an illustrative example, where 

x = 6 

■ ik x 10 5 6 

k 


v =0.3 


The results are plotted in figure 13 over the circumference of a cross 
section through the fuselage. 

This diagram shows the following features which are of practical 
interest: 

(1) The disturbance produced by the door opening is restricted to 
a rather small part of the shell. At an angular distance of 30° from 
the edge it has practically vanished. 

(2) The disturbance in the hoop force is without importance. It 

is only slightly higher locally than in the undisturbed part of the shell. 

(3) There are considerable stresses in the x-direction. The forces 
N x shown in the diagram are, of course, additional to forces which may 
exist from other causes. In particular, there is a zone of tensile stress 
near the edge. When taken together with a compressive force in the 
adjacent bar of the door frame, these stresses are comparable with bending 
stresses in a beam of span 2, which receives the load = pa from the 

undisturbed shell and is supported on the two rings shown in figure 12. 

(^) The forces N x are arranged in alternating tension and compres- 
sion zones of approximately equal width and decreasing intensity. The 
width is such that in usual fuselages it may be of the same order as the 
distance between stringers. If stringers were placed at the zeros of 

they would not influence our problem. This justifies the procedure used, 
which is based on the assumption of an isotropic shell without stiffeners. 

(5) If an additional stringer were provided right at the peak of N x , 

the stress distribution would, of course, be changed considerably. The 
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essential effect of this measure would be beneficial: The stringer would 

share the axial load with the shell. Therefore , it may be suggested that 
only the width of the first tension zone be determined, and a good stringer 
be provided at its center, that is, at a distance of 90 °/ K i from the 

door frame . 


Bulkheads 


General formulas for shells of revolution . - A cylindrical pressure 
cabin must be closed at its end by a bulkhead. This bulkhead may be 
constructed as a' flat wall built up of vertical and horizontal beams and 
a metal sheet. The beams have to transmit the air pressure by bending 
stresses' to the circumference of the bulkhead, from where it can be 
transferred to the cylindrical cabin wall. Since the total air pressure 
on the bulkhead is a force of considerable magnitude, a flat bulkhead 
will result in a heavy construction, which should be avoided if possible. 

The preferable shape of a bulkhead is that of a shell similar to a 
boiler end. When the cabin has a circular cross section, such a bulk- 
head will be a shell of revolution. As a basis for its stress analysis 
a short account of the theory of such shells will be given here. , 

Figure l4 shows the middle surface of a shell of revolution; its 
intersections with planes normal to its axis are parallel circles, and 
its intersections with planes containing the axis are all equal to each 
other and are called meridians. At all points of a parallel circle the 
angle 0 between its plane and a tangent to the meridian has the same 
value and is therefore characteristic for this circle. The angle between 
the plane of a meridian and the vertical will be called 0. Since a 
point of the shell is -determined by the parallel and the meridian on 
which it lies, the angles 0 and 0 may be used as coordinates on the 
shell. 

If the shell is cut along a parallel circle (fig. 15) , the stresses 
transmitted there can be found. As is usual in shell theory, equations 
are not written for the stress but for the meridional force N 0 which 

acts on the unit length of the circle. This force has the direction of 
the meridian. The resultant of all the meridional forces acting on one 
parallel circle is horizontal and of the magnitude N^(sin 0)2«r, and it 

must be equal to the resultant R = itr%> of the air pressure. Hence 



R 

2 «r sin 0 


pr 

2 sin 0 


(lfla) 
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When the shell is cut along a meridian, internal forces, that is, 
the hoop forces Ng, are found, hut they are not the same at all points 
of the- meridian and, therefore, cannot he found as simply as N 0 . A 

shell element limited hy two adjacent meridians and two adjacent parallel 
circles (fig- l 6 ) has to he cut out. The sides of this element, which 
are parts of meridians, have the length r^_ d 0 , where r-^ is the radius 

of curvature of the meridian. The other two sides have the length r d© 
(slightly different from each other because r is not the same on hoth 
parallel circles). The equilibrium of the forces N 0 (r d©), % ( r i 

and the air pressure p(r d©)(r-j_ d 0 ) in the direction of a normal to the 
shell yields the equation 


N 0 (r d©)d 0 + NgC r^ d 0 )d© sin 0 = p(r d©)(r-|_ d 0 ) ' 


The factor sin 0 in the second term comes from the fact that the 
resultant of the hoop forces lies in the plane of the parallel circle 
and has to he projected on the normal to the shell. The equation may 
he simplified to 



P 


Introduc ing N 0 


from equation (l 8 a) into this equation, 


% = 


Pr- 


sin 0 y 2 r 1 sin 0 y 


(18b) 


Equations (l 8 a) and (l 8 b) are sufficient to find the internal forces N 0 
and Ng when the shape of the shell is known. In order to permit the 

best use of the space in the pressure cabin and in the fuselage at its 
rear, an ideal bulkhead should he as flat as possible. This might lead 
to a bulkhead designed to meet the cylindrical cabin wall at an angle, 
as shown in figure 17(a). At the edge of the bulkhead equation (l 8 a) will 
yield a certain value of the force N 0 which, of course, must have the 

direction of a tangent to the meridian. Now this force cannot be trans- 
mitted to the cylinder because this shell can only resist a force 

having the direction of a generator. The difference, that is, the radial 
component of Nm, must be transmitted to a stiffening ring. It leads 
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there to a compressive force of considerable magnitude. The corresponding 
deformation, a decrease of the ring diameter, fits in no way to that of 
the cylinder and no better to that of the bulkhead. Therefore, all the 
trouble with bending stresses described in the section entitled "Bending 
of a Cylindrical Shell" arises here again but in a much more serious 
magnitude. The meridian of the bulkhead should, therefore, always end 
with a tangent parallel to the generators of the cylindrical cabin wall. 

At the center of the bulkhead r = 0 and 0=0, and formulas 18(a) 
and. 18(b) become indefinite. If the meridian has a finite radius of 
curvature r^ at this point, then, in its vicinity the relation 


r = r^ sin 0 

holds. Introducing it into equations (l8), they yield 


N 0 = % = ( l8c ) 

The tendency to make the bulkhead as flat as possible might lead to 
a meridian with an extremely feeble curvature in the central part. In 
the extreme case, for the curvature l/r-^ = 0, the stresses become infinite. 

This is illustrated by figure l8. The meridian in the upper half is a 
biquadratic parabola a 8 x = r , and the diagrams of the forces N0 
and N@ show the consequences of insufficient curvature of the shell. 

The lower half of the figure shows how easily the situation can be 
improved. Here the central part of the shell is replaced by a spherical 
segment, and at once the stresses are reduced to a moderate magnitude. 

It may be mentioned that the meridian chosen for this example does 
not fulfill the condition of a smooth transition to a cylinder, and, 
therefore, cannot be recommended even in the improved form, but the 
essential effect shown in figure 18 is, of course, true for any other 
shape with insufficient curvature. 

Ellipsoidal bulkhead .- An oblate ellipsoid (fig. 19) used as a 
bulkhead provides a good compromise between the desire to avoid dead 
angles and dark comers and the necessity of providing a smooth flow 
of forces. Relations between the radii r, r^_, and 0 may be found 

from the equation of the elliptic meridian. They are: 
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r 


sin 0 


(a 2 sin 2 0 + b 2 cos 2 0) 1/2 


r l 


(a 2 sin 2 0 + b 2 cos 2 0)^/ 2 


When these are introduced into equations (l8a) and (l8b), 

pa 2 

^ 2^a 2 sin 2 0 + b 2 cos 2 0j^/^ 

pa 2 b 2 - (a 2 - b^) sin^0 
2b 2 ^ a 2 s i n 20 + ^ 2 C os 2 0)^/^ 

These formulas describe completely the stresses in the bulkhead . They 
are not limited to a shell of constant thickness, hence the local 
stresses may be found simply by dividing by the local thickness t of 
the shell: 


a f = *0/* 

a e = ®e/ t 

The stresses at tvo points are of main interest: The center 0=0° 
and the edge 0 = 90°. In the center is found a biaxial tension 


°0 = a e 


2bt 

which determines the wall thickness. At the edge, the force trans- 
mitted to the cylindrical fuselage is independent of b and is the same 
for all ellipsoids (and for any other shape of bulkhead with smooth 
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transition to the cylinder). The hoop force N e depends largely on the 
ratio a/b of the axes. If b = 0.707a it becomes zero, and if the 
ellipsoid is still flatter the hoop force will be a compression. Since 
it is desirable to build the bulkhead as flat as possible, this fact 
deserves special attention. The compressive stress may be rather high, 
but it is confined to a small zone. Figures 20(a) and 20(b) show two 
examples of the stress distribution. In any case it will be wise to 
provide for a stiffening ring at the connection between the bulkhead and 
the fuselage. 

The compressive hoop stress has still another consequence which 
needs consideration. It produces an elastic deformation, which decreases 
the diameter of the boundary circle of the bulkhead. On the other hand, 
the diameter of the cylindrical wall of the fuselage will increase in 
the part in front of the bulkhead as a consequence of the positive hoop 
stress in a cylindrical shell, and will not change at all in the part 
behind the bulkhead where there is no internal pressure. The deforma- 
tions of the three shells look somewhat as shown in figure 21. Since 
the shells are connected to each other, such a discrepancy cannot exist 
in reality but will be prevented by a system of bending stresses in the 
boundary zones of all three parts. In boilers and other pressure vessels 
these bending stresses are rather serious and means to avoid them are 
desirable. In pressure cabins they may be of some minor importance, but 
It is certainly better to eliminate them as far as possible. 

The discrepancy between the two cylindrical parts is, of course, 
unavoidable, but the edge deformation of the bulkhead should lie between 
those of the two cylinders. That means, at least, that the hoop force 
must be positive. A hemisphere would fulfill this condition per- 
fectly, but as a bulkhead it would lead to poor utilization of space. 

A better solution is to turn the ellipsoidal bulkhead with its convex 
side toward the pressure cabin (fig. 22(a)). Then the previous for- 
mulas are still applicable, but the signs of all stresses are reversed. 
That is desirable at the edge but certainly not in the center, where 
compressive stresses create a buckling problem. This will be avoided 
by the bulkhead shown in figure 22(b). The trouble with this shape is 
that the membrane stresses at the sharp edge between the convex and 
concave shells cannot make equilibrium with each other without the help 
of a stiffening ring. This is shown by figure 23. The two meridional 
forces N0 shown there, one a tension and the other a compression, will 
have the same horizontal component and thus assure the axial equilibrium 
of the shell; but their radial components are both directed outward and, 
therefore, cannot equilibrate each other. If a ring is provided, they 
may both be transmitted to it as shown, leading to a compressible hoop 
force in the ring. But here a new difficulty arises. The hoop strains 
of the three parts will not fit together and will again lead to bending 
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stresses. It may be that they will turn out to be less serious in a 
particular case, but at least they are now at a place where they do not 
produce quilting of the surface of the fuselage. 

Bulkheads with improved boundary effects .- A more promising method 
to avoid excessive bending stresses would be to choose another shape of 
the meridian. The fact that a hemisphere gives just what is wanted, an 
edge deformation halfway between zero and that of the cylindrical cabin 
wall, indicates that the meridian must begin with a curvature l/a at 
the edge. In order to make the bulkhead flat, the curvature should then 
increase and later become very small when the center is approached. 

This is schematically shown in figure 24, but the idea cannot be executed 
in this form because any discontinuity in the curvature again will pro- 
duce those local bending stresses that are to be eliminated. 

What is needed is a curve having the same general shape but a smooth 
transition of curvature. Such a curve may be found in several ways. 

One is a modification of the Cassinian curves (see also reference 5). 
Using the coordinates x and r in a meridional plane as shown in fig- 
ure 24, the equation 

(n 2 x 2 + r 2 ) 2 - 2A 2 (n 2 x 2 - r 2 ) = B^ ( 19) 

describes a Cassinian curve, if n = 1. With n > 1 the curves are 
flattened and for 2 < n < 3 assume a reasonable shape. The parameters 
A and B must be chosen so that for x = 0 the ordinate becomes r = a 
and the radius of curvature is r^ = a. This yields 


A = 


a 



- 1 
+ 1 



The radius of curvature at the center of the bulkhead will then be 


r 


1 




- 1 

+ 1 


a 
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With these data the internal forces at the most interesting points 
can he found. At the edge equations (l8a) and (l8h) with r = r^ = a 

and 0 = 90° give: 




pa 

2 


At the center equation (l8c) must he used and 


N0 = N e 



1 



- 1 
+ 1 


is obtained. 

But this is not enough. To he safe from surprises, one must have 
the stress 'distribution along the meridian. It can he found by the 
following procedure. 

* 

Assume x anH from equation (19) find r, or vice versa, depending 
on which will give the greater accuracy . Then compute 


pop p p 

n x n x + v - A^ 
r n 2 x 2 + r 2 + A 2 

— = r ^1 + (r') 2 
sin 0 

*i d 2 r 

r -J? 

r' , (r* n 2 x 4A 2 (n^x + rr 1 ) 

x r r (n 2 x 2 + r 2 + A 2 ) 2 
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r^ sin 0 1 + (r 1 ) 

Nov formulas (l8a) and (l8b) may be applied. 

This was done for n = 2 and the results shown in figure 25 were 
obtained. The hoop force Ng falls to zero and rises at the center of 

the bulkhead toward the value previously mentioned. If n > 2 is 
chosen, the hoop force will become negative, and very much so if n is 
too great. This, of course, should be avoided. 

Bulkhead for double cylinders . - In a double- cylinder cabin the two 
cylinders may have their bulk h eads at different stations. Between the 
two bulkheads the longer cylinder must have a full circular cross sec- 
tion, and its' intersection with the first bulkhead leads to a difficult 
stress problem. 

This situation is eased considerably if it is possible to have both 
bulkheads at the same cross section. Their shapes may then be chosen 
such that they intersect in a plane horizontal curve. This will be the 
case when they are oblate ellipsoids with the ratio bra = b’ra' . They 
may then be derived by affine transformation from two spheres as indi- 
cated by figure 26. 

A reinforcing ring must be provided along the intersection of the 
two shells. It will now be shown that this ring, which has the shape 
of a half-ellipse, will be stressed in its own plane Only, and, with a 
certain exception, will even be free from bending stresses. 

In the case of two spheres this is evident. Introduce coordinates 
0 , 0 ' , and 6 as shown in figure 26. If the two shells are inflated 

by a pressure p, the internal forces will be 

N0 - | pa 

in the upper shell and 

y - v - 1 

in the lower. The forces for 0 = l80° - a and N^' for 

0 ’ = l80° - a' act on the ring as shown in figure 27. Since 
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a sin a = a' sin a' 

= c 

the vertical components N 0 sin a and. N^' sin a' balance each other 

and the horizontal components combine to a uniform radial load, pro- 
ducing a positive hoop force 


F = c i p(a cos a + a 1 cos a' ) 


in the ring. 

When there are two ellipsoids, the stress analysis is not so simple. 
One might think of using the solution given under the section entitled 
"Ellipsoidal Bulkhead" for a single ellipsoid and of determining from it 
the forces acting on the elliptic ring. But there is no reason to believe 
that the axial symmetry assumed in deriving that solution still exists 
when part of the shell has been cut away and its symmetry thus destroyed. 
The clue to the solution is the idea that the ring shall be free of 
bending moments, and it has only to be shown that a solution with this 
property exists, that it is unique, and how to find it. 

To fulfill this program, some notions of the theory of affine shells 
are needed (for details see references 2, 3 , and 5). They will be pre- 
sented here as applied to the particular problem to be solved. 

In addition to the curvilinear coordinates 0,0 and 0‘,0 on the 
two spheres, two systems of rectangular coordinates are now introduced: x, 

y, and z for the ellipsoids and x*, y*, and z* for the spheres. The 

simple geometrical relation between the two shells is represented by 


x = nx* 



with n = b/a. As curvilinear coordinates on the ellipsoids the 
values of 0 (or 0 *) and 0 for the corresponding points on the 
spheres are used. These coordinates do not represent angles that can 
be measured on the ellipsoids, but each pair of values 0,0 defines 
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clearly one point on the shell, and. this is all that coordinates are 
expected to do. 

A shell element is cut out of the upper sphere by two meridians 6 
and 6 + d0 and by two parallels 0 and 0 + d0. It has the area 

dA* = a 2 sin 0 d0 d 6 

When it is projected on the planes (y*,z*), (z*,x*), and (x*,y*), the 

projected areas are found: 

cLA x * = dA* sin 0 cos 6 

dA y * = dA* sin 0 sin 6 

dA z * = dA* cos 0 

The element of the bulkhead shell is simply the projection of the 
element dA* on the ellipsoid. Both have the same projection on the 
yz-plane,' but the other two projections are reduced in the ratio n:l: 

dA x = dA x * 

= dA* sin 0 cos 6 

dAy = n dAy* 

= ^ dA* sin 0 sin B 
a 

dA z = n dA z * 

= — dA* cos 0 
a 

When the ellipsoid is inflated by an Internal pressure p, the force 
acting on the element is p dA and has the rectangular components p cLA x , 
p dAy, and p dA z parallel to the axes x, y, and z. This load pro- 
duces internal forces transmitted at the four sides of the shell 
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element dA. On each side of the element this force lies in a tangential 
plane to the shell and may be resolved into a normal and a shear com- 
ponent. Something better can be done: Oblique components parallel to 

the coordinate lines 0 = Constant and 9 = Constant can be used. 

These forces divided by the length ds^ or dsg of the line element 

are called N0, Nq, and N as indicated in figure 28. 

A simple relation between these forces and a certain system of 
internal forces N0*, Nq*, and N0 e * in the sphere will now be 

established. These forces must, of course, lie on tangents to the 
sphere, and they will be chosen in such a way that they have the same 
projections on the yz-plane as the corresponding forces in the ellipsoid. 
Then both will have the same components in t the directions y and z, 
but the x-components of the forces N in the ellipsoid will equal n 
times the x-components of the forces N* in the sphere. 

The forces N* of the sphere will be in equilibrium with a load 
which has the same y- and z-components as that applied to the ellipsoid, 
but l/n times its x-component. 

If the load components per unit area dA* of the sphere are denoted 
by P X *> Py*> and P Z *> 



= p £ dA* sin 0 cos 6 

D 

p y * dA* = p dAy 

= p dA* sin 0 sin 6 


p z * dA* = p dA z 

= p ^ dA* cos 0 


/ 



34 


NACA TN 2612 


and hence the loads per unit area of the surface of the sphere are 


p * = p — sin 0 cos 6 
■* b 

Py* = p ^ sin 0 sin 6 
Pz* = P | cos 0 


( 20 ) 


If the forces set up in the •Spherical shell by this load can be found, 
it is only necessary to project them on the tangential plane of the 
ellipsoid and to refer to the unit length of the ellipsoid's line element 
and then the forces in the ellipsoidal bulkhead vill be obtained. 


To make the stress analysis for the sphere along conventional lines, 
the loads given by equation (20) are transformed into the components 
X*, Y*, and Z* as shown in figure 28: 


X* = -p * sin 6 + p * cos 6 
A y 


p (a 2 - b 2 ) , 

— £in 0 sin 26 

2ab 


Y* = ^P x * cos B + Py* sin e) cos 0 - p z * sin 0 
p(a 2 - b 2 ) 


2ab 


cos 0 sin 0(l + cos 26) 
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These expressions have the form 


X* = X Q * + Xg* sin 26 
Y* = Y 0 * + Y 2 * cos 26 
Z* = Z Q * + Zg* cos 26 

showing the harmonic constituents of orders 0 and 2. 

The first term on the right is a load with rotational symmetry: 

V = 0 

! 

p(a 2 - b 2 ) . A 

V = cos 0 sin 0 


Z o* = 2ab + ^ ~ sin \ 


Simple formulas , which will not he reproduced here, lead to the internal 
forces: 


Y - N0 O * 

- £*> 

2 


V - n eo * 

= |^(h 2 + ( &2 - h 2 )sin 2 0j 


( 21 ) 


N 0e* = 0 
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These formulas are valid for the upper sphere of radius a. For the 
lower sphere it is necessary to write simply a', b', and 0' , instead 

of a, b, and 0, respectively. 

The second harmonic of the load. 



p(a2 - b 2 ) 
2ab 


sin 0 


p (a 2 - b 2 ) 


2ab 


cos 0 sin 0 


* p (a 2 - b 2 ) o , 

Zo = — ' sin 2 0 

2ab 

may not be handled so easily. It leads to forces which depend also on 
a sine or cosine of 28 and may be written as 


N 0* = N 02* cos 26 


V = N 02* C0S 26 


N 0e* = N 062* sin 26 

The basic formulas connecting N^* , ^62*^ an ^ ^062* with the load- 

components % 2 * > ^ 2 *) and Z 2 * may be found in the literature (refer- 
ence 2, pp. 37 to 44, or reference 3). They lead to a solution having 
two free constants for each of the spheres. One constant in each pair 
may be determined from the condition that the stresses are finite at 
0=0 and 0' =0. The other two are still to be determined. In this 
way are obtained 
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W 02 


* p(a c - b 2 ) 


213 


2(1 + cos 0 )‘ 


N 002 


* p(a 2 - b 2 ) ^ C 

- — i cos 0 - 

213 2 (l + cos 0) 2 


N, 


02 ~ 


1 (a 2 - b 2 ) 


2 b 


cos 2 0 


2(1 + cos 0 )' 


for the upper shell and corresponding formulas for the lower, containing 
another constant c' . 


From these forces the ring receives a radial load 
(N 02 * cos a + ^ 02 *' cos a ') cos 2 ®> a vertical load 

(N 02 * sin a- - n 02 *' sin a ') cos 2e > and - a shear l° ac i ( N 002 * + N 002 *') sin 26 •» 
positive as shown in figure 29. The two free constants give the opportunity 
of influencing these forces in such a way that the ring is free of bending. 
The first thing to do is to make the vertical load vanish at every point. 
This yields the equation 

N 02 * sin a - N 0 2 *' sin a* = 0 

When the express ion just given for N0 2 * is introduced and some simple 
geometric relations mentioned before are used, the following equation 
is obtained: 



sin a 
- cos a ) 2 


C' 


sxn a' 


(l - cos a ' ) 2 


( 22 a) 


The second equation in C and C' must express the fact that there 
is no bending in the plane of the ring. Under this condition the ring 
has only an axial force F, which, of course, will be a function of 0. 
From the equilibrium of the ring element (fig. 30 ) two relations are 
found: 


c(l02* cos a + n 02*' 


cos a 


') 


cos 20 = F 



N002*' ) sin 26 = 


dF 

d0 
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From the first one it is possible to write F = F£ cos 26, and, elimi- 
nating F£ from both equations, 

N0 2 * cos a + w 02*' cos a * + !( N 082* + N 082*') = 0 

When the expressions found for the internal forces of both shells are 
again introduced and then simplified by appropriate use of geometric 
relations, this will yield the second equation for C and C*: 

(1-2 cos a) , (l - 2 cos a') 3p(a 2 - b 2 ) sin (a + a' ) , 

C — + C' — = ; (22b) 

(1 - cos a) (1 - cos a’ 

These two equations, when solved in general terms, yield 


C = 


3p(a 2 - b 2 ) (l - cos a) 2 sin (a + a 1 ) 

b sin a + sin a* - 2 sin (a + a' ) 


Introducing this into the formulas for the internal forces, the following 
expressions for the upper sphere are found: 


N 02 


p(a 2 -b 2 ) 


2b 


1 + 


3 sin (a + a' ) 


sin a + sin a' - 2 sin (a + a' ) 


(l - cos a) 2 
(l + cos 0) 2 


N 092* = 


p(a 2 - b 2 ) , 

i + 

3 sin 

(a + a' ) 

— 

2 

(l - cos a) 

<c.uo y 
2b 1 

U T 

sin a + sin a' 

- 2 sin 

(a + a' ) 

(l + cos 0) 2 


y (23) 


p(a 2 - b 2 ) | 

J. 

3 sin (a + a’ ) 

(l - cos a) 2 

a. 1 

CUo yj *r 

sin a + sin a' - 2 sin (a + a' ) 

(1 + cos 0) 2 


These and those in equations (21) are the forces set up in the 
spherical shell by the fictitious load (20). The last step necessary 
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is to find from these forces the real forces in the ellipsoidal bulk- 
head under the uniform air pressure p. 

The center of the bulkhead has the coordinates 0 = 90 ° , 0=0°. 

The internal forces are here parallel to the yz-plane and, therefore, 
are the same in the sphere and the ellipsoid: 


N 0 - N 0O* + N 02* C0S ° C 


pb p(a 2 - b 2 ) 


2b 


1 + 


{3 sin (a + a*)J(l - cos a) 2 
sin a + sin a’ - 2 sin (a + a*) 


pa. 

2b 


2 , 3p(a 2 - b 2 ) [sin (a + a*)](l - cos a)‘ 


2b sin a + sin a* - 2 sin (a + a’) 


N 0 = N 0O* + N 02* cos °° 

pa 2 3p(a2 _ b 2 ) [sin (a + a*)J(l - cos a) 2 
2b 2b sin a + sin a* - 2 sin (a + a*) 


N 00 ~ N 00O* + N 002* sin °° 

= 0 

■yt- 

* 

On the edge 0 = ± 90 °. Here the force Wg Has x-direction and 
must be reduced by multiplying by n = b/a. The force N 0 is parallel 
to the yz-plane; therefore, the forces on corresponding line elements 
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of sphere and ellipsoid are the same hut the line element is reduced 
in the ratio b/a, hence the force per unit length of this element 
increased by a factor a/b. The shear is zero. Thus on the edge of 
the ellipsoid: 


>V ■ ?(.V + >2* C0S l8 °°) 


3 sin (a+'a*) (l-cos a) 2 

” ‘ 2 
sin a + sin a* - 2 sin (a + a* ) (l + cos 0) 

p(2b 2 - a 2 ) a 3pa(a 2 - b 2 ) __ sin (a + q 1 ) (l-cos a) 2 

2b 2 2b 2 sin a + sin a* - 2 sin (a+ a 1 ) (l+ cos 0) 2 


pa pa (a 2 - b 2 ) 
2 2b 2 


w e = |( N eo* + %2* cos l8o °) 


pb 2 p(a 2 - b 2 ) 


2a 


2a 


O . , 

sin 0 + cos^0 + 


3 sin (a + a* ) (1 - cos a) 


sin a + sin a’ - 2 sin (a + a’ ) ^ + cos qy 


P& 3p(a 2 - b 2 ) sin (a + a 1 ) (l - cos a) 2 

2 2a sin a + sin a' - 2 sin (a + a' ) ^ + cos 


An example for the distribution of these forces over the circumfer- 
ence is shown in figure 31. In the top and bottom zones the forces are 
rather uniformly distributed, but there are marked peaks at the junction 
of the two ellipsoids. This fact indicates that stresses in shells 
should be determined carefully. 
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In all the preceding formulas, the denominator 


sin a + sin a' - 2 sin (a + a* ) 


appears. It may happen that a and a.' have such values that this 
denominator is zero, in which case the formulas would yield infinite 
stresses. This indicates that in such a case the combination of the 
two spherical or ellipsoidal shells is capable of an inextens ional 
deformation and that rigidity can only be secured by giving the neck 
ring sufficient rigidity against bending in its own plane. If plane 
bending of the ring is to be assured, then equation (22a) between the 
two constants C and C’ still must hold. The bending moment In the 
ring then becomes independent of the choice of C. Equation (22b) 
becomes useless and must be replaced by the condition that the internal 
forces assume finite values. This leads to C = C ! = 0. 

The somewhat lengthy analysis of the double bulkhead has been 
reproduced here not only because of the particular problem under con- 
sideration but as an example of two important features of thin shells: 

(1) The fact that the stiffening ring along the intersection of 
two parts of the shell is usually free of bending moments 

(2) The use of affine relations for the solution of shell problems 


Nose of Plane 

General rules .- The nose of the fuselage may have so many various 
shapes that not much can be said in general about its stress analysis. 

In high-speed planes aerodynamic consideration may lead to shaping the 
nose as a perfect surface of revolution. If it is part of the pres- 
surized cabin, it may be treated with formulas (l8a) and (l8b) for 
stresses in such she 11s ; if the cabin terminates in a bulkhead back of 
the nose, all that has been said about the rear bulkhead is applicable. 

The modern passenger plane usually has a nose which looks like that 
shown in figure 32. The major part of it is a shell, but the smooth 
surface is interrupted by many windows. In such cases a shell analysis 
as described in the preceding sections will, in general, be too com- 
plicated for practical purposes. As regards the stress analysis of such 
structures, the following facts should be kept in mind: 

(l) All large uninterrupted parts of .the metal skin will act as 
shells, whether they are fixed on a solid framework or only stiffened 
by rings and stringers. The stiffeners which are connected to the shell. 
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although absolutely necessary for the introduction of local loads and 
as a buckling reinforcement, are obstructions to a smooth flow of stress 
in the shell proper and lead to quilting and to tensile stresses in 
rivets. 

(2) All edges of such shell parts, for example, along the windows, 

must be stiffened by edge members. It is always advantageous to shape 
these edge members after plane curves. With rare exceptions they will 
not then be subjected to bending in space, but they must offer resistance 
to bending in their plane and require the corresponding rigidity, bracing 
and support. . • ' 

(3) There is no need for making cross sections circular. Any curved 
shell can resist an internal, pressure, but, of course, the stress dis- 
tribution will be less uniform and may easily have local zones of com- 
pression if the cross sections are far. from circular. 

(4) Areas of extremely low curvature should be avoided. Membrane 
shell theory leads to extremely high stresses in such parts and, owing 
to these stresses, the panel bulges out, thus increasing the curvature 
and reducing the stress. In addition this bulging invariably leads to 
some plastic- Jieformation at the edges of the panel and, therefore, to 
a permanent bulging, -which is undesirable. 

Ellipsoid with three different axes .- A general ellipsoid is a kind 
of a shell which probably will not occur as part of a pressure cabin.- . 
However, its membrane forces are easily computed and may give an idea 
of what may happen in other shells of noncircular cross section. 

Consider an ellipsoid having the three half-axes a > b > c and 
being subjected to an internal pressure p. In order to find the 
membrane forces, establish relations between them and those in a sphere 
of radius b under a certain load. This follows the same lines as the 
theory for the double bu lkh ead in the section entitled "Bulkhead for 
Double Cylinders . " . _ 

In rectangular coordinates x*, y*, and z* the sphere has the 

equation . 

x* 2 + y* 2 + z*^ = b 2 

and in coordinates x, y, and z the equation of the ellipsoid is 
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As surface coordinates on the sphere the angles 0 and 0 are used 
as shown in figure 33- Through the relations 



each point of the' .sphere corresponds to a point on the ellipsoid. By 
attributing the same values of 0 and 0 to both, a system of coordi- 
nates is established on the ellipsoid. Its lines 0 = Constant are 
parallel ellipses in horizontal planes. Its lines 0 = Constant are 
ellipses in planes through the z-axis. 

The shell element on the sphere has the area 


dA* = a d0a sin 0 d0 

and its projections on the coordinate planes are 

dA x * = dA* sin 0 ' cos 0 
dAy* = dA* sin 0 sin 0 
dA z * = dA* cos 0 

The projections of the corresponding element of the ellipsoid are 

- dA x = § dA x * 

• dAy dAy* 

. - 4 -dA 8 -f'dA 2 * 
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Multiplying this by p yields the components of the force p dA acting 
on that element. The corresponding forces on the spherical element are 
then 

p x * dA* = |-p dA x 

= ^ p ^ dA* sin 0 cos 8 

Py* dA* = p dAy 

= p dA* sin 0 sin 8 
p 2 * dA* = | p dA z 

' = ^ p ^ dA* cos 0 

From these relations are found p x *, p y *, and p z *, the loads per unit 

area of the sphere, in directions x*, y*, and z*. They are connected 

with the usual components, X* in direction 8 , Y* in direction 0, 
and Z* in the radial direction by the formulas: 

X* = -P x * sin 6 + Py* cos 8 

Y * = (p x * cos 8 + Py* sin e) cos 0 - p z * sin 0 

Z* = ^p x * cos 8 + Py* sin &) sin 0 + p z * cos 0 
Introducing p x *, p y *, and p z *, they yield: 

** ■ I £(£ - I) sln * 5ln 26 

= X 2 sin 28 
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Y* = ^ cos 0 sin 0 


r £ , ac a 
+ h 2 " * c 


c/a b\ 
b \b " a / 


cos 26 


= Y q + Y 2 cos 26 



2 a + 

c 



ac 

b2 




= Zq + Z 2 cos 26 


The corresponding membrane forces may be found from well-known formulas 
(reference 2, pp. 37 to 39). They are: 



To find the forces in the ellipsoid, the ratio of the corresponding line 
elements is needed. The results are mentioned only for the three points 
A, B, and C. 

At the point A 0 = 90°, 6 = 0°, and 
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to - P^/W . ac al> \ 

0 2c \ a b c / 


p /b^ ab £ 


= Pa /, h£ 

2 c2 a 2 ; 


At the point B 0 = 90°; © = 90°, and 


pb 


C 2 c 2 


N 0 " 2 l 1 + b 2 a 2 


% 





At the point C 0 = 0° and 6 may have any .value. When 0 = 0° is 
assumed arbitrarily, then lies in the xz-plane, and Ng, in the 

yz-plane, and the internal forces are: 




If (a2 + b^)c^ < a^b2, the force % at the points A and B becomes 

a compression. The other four formulas always yield positive forces 
when a 2 b ^ c. These results may serve as a first orientation of what 
is to be expected in shells of noncircular cross section under the 
action of an internal pressure. 



NACA TN 2612 


47 


PLATE PROBLEMS 
Stresses in Thin Flat Sheets 


Plane plates are not very desirable as parts of the wall of a pres- 
sure vessel, but it often is not possible to avoid them. Therefore, the 
stresses set up in them by a lateral pressure p will be considered 
here. 


If such a plate were thick enough, it might carry its load by 
bending stresses as does the reinforced concrete floor slab of a 
building; however, the skin of an airplane is much too thin to carry 
an appreciable load with tolerable bending stresses. Its stress system 
is a superposition of bending stresses and of the stresses in a flexible 
skin. 


The subject of this section will be such a thin skin of rectangular 
shape. Its stress problem is essentially nonlinear. In two dimensions 
it is so involved that all theoretical and experimental effort spent 
on it up to the present is still far from giving a complete answer to 
all questions which the engineer might ask. Therefore, a discussion is 
first presented for the one -dimensional problem which, in many cases, 
will give useful information for practical purposes and beyond that will 
show the general features of the stress system present in the two- 
dimensional case. 

Thin sheet stressed in one dimension .- Consider a thin plate as 
shown in figure 3^. In the x-direction it has the span Z, and the 
'Sides x = tl/2. are supported in such a way that not only the deflec- 
tion w but also a displacement u in the x-direction is prohibited. 

In the direction of the y-axis the plate is supposed to be long enough 
to make the conditions at the shorter sides immaterial. 

For the purpose of stress analysis cut a strip of unit width out 
of this plate. Because of the end conditions, the lateral load p will 
produce a direct stress a along the strip, which is necessarily 
independent of x. If the deflection w is large enough, this stress 
will be capable of carrying the load. 

This can be seen "on an element of length dx cut from the strip 
(fig.. 35). The condition of vertical equilibrium is: 

/dw\ 

at df 1 + p dx = 0 
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It yields the relation 

d^w _ p 
dx 2 crt 

which indicates that the strip must deflect into a common parabola: 


p* 2 L ^ g \ 

8crt\ “ z 2 ) 


•With a maximum at x = 0: 



(24) 


The horizontal displacement in the direction of increasing x may 
be called u. The strain in the strip is then 


e 


du + 1 /dw\ 2 
dx 2 \dx/ 


a 

E 


At the center x = 0, u = 0 from symmetry. At the support 
x = Z/2, therefore, 



nZ p 2 Z 3 
33 " 48a2t 2 


(25) 


Since an unyielding support u Q = 0 was assumed, this equation may be 
used to find a for a given load: 


l v 3 Ep 2 Z 2 

2 \|~ 


a 
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From equation (24) the final expression for the deflection may he found, 
in particular that for the maximum w = f at x = 0: 


l ,3f3pT 

4 VEt 


These two formulas represent nonlinear relations, owing to the quadratic 
term which represents the influence of the deflection w on the strain e. 
If this term should he neglected, as it is in many other cases, no 
reasonable result at all would he reached. The nonlinearity is, there- 
fore, an essential feature of this problem. 

Since the plate deflects, there will he a bending moment 





? t3 P 

E ii « 



This is incompatible with the assumed support and has been neglected in 
the preceding formulas, but this may safely be done if M is small as 
compared with the moment pZ 2 /8 which would be necessary to carry the 
load by beam action. This condition may be brought into the dimension- 
less form: 


For a pressure cabin with p 
this yields 



_3_ 

16 


= 7 psi and duralumin with E = J.0^ psi 
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For the airplane, the assumption of unyielding supports goes too 
far. The edges x = Constant of the plate are kept apart by stiffeners 
in the x-direction. When they are not riveted to the skin, the problem 
is still one -dimensional and may be represented by a flexible strip and 
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a strut as shown in figure 36. The flexible strip is exposed to the 
load p and deflects under it; its ends are kept apart by a strut of 
cross section A]_, which has a compressive force K = -crt; the supports 

are such that they allow the corresponding elastic deformation of the 
bar. 


When the cross section of one stiffener is A , and the distance 
between stiffeners is d, then ihe area Aj_ = A/d corresponds to a 1 

strip of unit width of the plate. 

When the horizontal displacement u is assumed to be zero at x = 
(midspan), equation (25) may again be used for the displacement u Q at 
the end, but now Uq must correspond to the fact that the strut becomes 
shorter by N2/EAj_: 

crtZ 

u 0 = - 

0 2EA 1 

\ 

Introducing this into equation (25), 


a 


1 y3 Ep^Z 2 .3 Ap 

2 \|aT^ 


(26a) 


for the stress in the plate. The greatest deflection follows then from 
equation (24): 


f = 


l 

4 



(26b) 


Comparing these two formulas with those which were obtained for 
nonyielding supports, it is seen that they become identical for A-^ — *■ 

Since t/Aj is more likely to be equal to 1 than 0, the assumption of 
nonyielding supports may lead to errors of about 25 percent, overrating 
the stress and underrating the deflection. It seems, therefore, not 
worth while to spend much effort on the two-dimensional problem if this 
effect is not taken into account. However, the simpler formulas are 
good enough for estimating the order of magnitude of a and f and 
for discussing the influence of the bending stiffness of the plate. 

The formulas (26a) and (26b) are sufficient if the sheet panel is 
part of a flat bulkhead. But in most other cases the wall, consisting 
of' the sheet and its stiffeners, has to transmit an internal force such 
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as N x or N 0 explained in the section entitled "Circular Cylinder. " 

These forces may he due to the over- all bending of the fuselage, to the 
action of the internal pressure on other parts of its wall, and to other 
causes. When the sheet is flat and bulges out, the distribution of this 
force on sheet and stiffeners is no longer governed by the formulas (5). 
It may be found by adding an axial force P to the strip and strut 
system of figure 36 (see fig. 37)* 

If a again indicates the stress in the sheet, the force in the 
strut is 


N = P - crt 


and the horizontal displacement at x = Z/2 must be 

U ° “ 2EA-L 

Equating this to the expression (25), 



_2 P_ _ Ep £l 2 
A 1 24t 2 


= 0 


(27) 


This may easily be solved in any given case. The deflection f follows 
then from equation (2b). 

Some results are shown in figure 38 in dimensionless variables. The 
values for the parameter pZ/Et have been chosen as rather extreme in 
order to cover the whole field of practical interest. For most of the 
diagrams, t/A]_ = 1 has been assumed, but one of them shows the trend 
for a variation of this parameter. ' 

The diagrams show the influence of the force P. They emphasize 
that a solution of the two-dimensional problem which disregards this 
influence cannot yield more than a rough approximation of the real air- 
plane problem, even if it were an exact solution of the simplified 
problem. 

Thin sheet stressed in two dimensions .- Assuming a sufficiently 
thin plate, the formulas developed in the preceding section are exact 
for infinitely long rectangles. It is probable that they will yield 
good results if the ratio of the sides is 1:4 or even 1:3, but when 
the rectangle approaches a square, they become inapplicable. 
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To find out how the results must he modified in such cases, consider 
a square plate framed by four equal stiffeners (fig. 39). When the plate 
bulges out under a lateral load p, stresses a will be set up in two 
directions, which are denoted ct x and Gy. The stiffeners will receive 

compressive forces. The plate in their immediate vicinity must have the 
same strain and hence a compressive stress. The distribution of Oy 

along one of the edges must, therefore, be as shown on the figure. 

At the center of the plate ct x = Gy. When two strips are cut out 

along the coordinate axes, each one will carry one -half of the load p, 
and the curvature will be half of that which would follow with this 
same a from the one -dimensional theory. When the horizontal strip is 
followed toward the right edge of the plate, the stress Gy will 

decrease and finally become negative. Where it passes zero, the one- 
dimensional theory will yield the correct curvature, and closer to the 
edge the curvature of the square plate will be greater than that of a 
single strip. A single strip was seen to deflect as a parabola. The 
profile of the square plate must, therefore, be far from a sine curve, 
and results computed on this assumption must, therefore, be interpreted 
with some reserve. 

When a diagonal of. the plate is followed a x always equals Oy , 

but both stresses decrease the father away they are from the center. 

The curvature in both directions must, therefore, become greater, and 
a sharp fold may be expected toward the corner. But at least there is 
a region where a x and Gy become negative and are no longer capable 

of carrying any load at all. Here even the thinnest plate must have 
essential bending stress. The thinner the plate is, the smaller this 
region will be, and the sharper the curvature will become. It follows 
that the highest stresses will occur on or near the corners. They may or 
may not be responsible for the ultimate load of the plate, depending on 
the possibility of smoothening the peak; by local plastic flow, but they 
are certainly responsible for permanent deformations which produce that 
quilting of flat panels which makes airplanes unsightly and is not much 
appreciated by the aerodynamic ist. 

There is little numerical information available on square and rec- 
tangular plates. Some papers (references 6 and 10) are mentioned in the 
references of this report. One of them, reference 6, contains rather 
complete material for plates with unyielding supports (fig. 34); however, 
this paper is based on the assumption that the profiles of the deflected 
plate along both axes are sine curves. In particular, this assumption 
is also made for the lengthwise profile of long plates (1:4), where it 
leads to an overrating of the influence of the support at the short 
sides. The results for the long plate, therefore, do not check with 
the one -dimensional theory. However, the two-dimensional problem is 
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so complex that a critical use of the diagrams of Moness (reference 6) 
is the test that can he recommended at this time. 


Thermal Stresses in Window Panes 

The window panes of a pressure cabin are not only exposed to the 
difference in pressure between the interior of the cabin and the free 
atmosphere but also to a considerable difference in temperature. The 
bending stresses due to the pressure may easily be found from text- 
book formulas, but the thermal stresses require some discussion. 

Consider a plate of uniform thickness, sjmply supported along its 
edge . . Assume that no load is applied but that there is a difference T 
between the temperatures of its faces. When the temperature is increased 
by T, a positive strain 


e = aT 


in every direction. 
When only one side 
strain between both 


in every direction. The middle surface of the plate is then deformed 
into a small part of a sphere of radius 1 / k . 

If the plate is circular, this is all that happens. The window 
will slightly deflect to the warmer side, and no thermal stresses will 
be set up. But if the plate is rectangular, the deformed shape will 
no longer fit on the support, and, in order to make it fit, the support 
will exert forces on the plate and these forces will produce bending 
stresses. 

Formulas for the bending moments which correspond to these stresses 
will now be established and discussed. In doing so, the following nota- 
tions will be used (fig. 1+0): 

x,y coordinates ' 

w deflection 

M x ,My bending moments 

M X y 


a being the coefficient of thermal 
of the plate is heated, e is the dif- 
sides, and this difference leads to a 


will occur 
expansion, 
ference in 
curvature 


twisting moment 
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K 


curvature in x-direction 



■bending stiffness 


Et 3 


U2 (l - v 2 ) ■ 


The stresses and deformation due to heating of the upper face of 
the plate will he built up in three steps. The first one has already 
been done, namely, the application of the temperature difference to the 
free plate, resulting in a uniform bending without stress. 

In the second step this deformation is completely removed by applying 
along all four edges of the plate, constant external bending moments M Q 

of appropriate size. They produce bending moments 


^ = My = Mo 

in the plate which are constant everywhere and in all directions. Now, 
the curvature k x of the plate is related to the bending moment by the 

well-known formula: 


= - v 

In this case, it yields 

s - - -Ha— 

K(1 + v) 

To remove the thermal deformation, « x must be made equal to -aT/t 
and therefore ' 


Mo = K(1 + V) 


Under the combined action of the temperature T and the edge 
moments Mq the plate is perfectly plane and may be attached to its 

supports. Now the third step may be done: The condition of simply 

supported edges must be realized and the edge load M Q compensated by 
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adding an edge load -M 0 . This is a problem of plate theory and may be 
solved in the following way: 

When there is no lateral load, the deflection w of the plate must 
satisfy the well-known differential equation (reference ll): 

S^w d^w S^w 

r- + 2 ■ ■ ■ + j- = 0 

8x^ cix^dy^ dy^ 

Introducing the sum of the two bending moments 


M = M x + My 


as an auxiliary variable, this equation may be split into two equations of 
the second order: 


8x^ 8y^ 


= 0 


(28a) 


c)^w _ M 

^2 + ^2 " ' K( 1 + v) 


(28b) 


These equations can be solved one after the other because a boundary 
condition can be found for each one. 

Consider, for example, the edge x = a/2 of the rectangle. The 
bending moment -Mq is applied there and 


- -«o 


is obtained. Since the edge is supported, w = 0 for all values of y 
and, hence, also 
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Introducing this into the elastic lav 




_ w 8^v\ 

^ ‘ w 5 v S * 2 / 


(29) 


it can be concluded that My = -VM Q on this edge and, hence, that 


M = -(1 + v)M c 


(30) 


A similar reasoning may be made for the three other edges. 

The solution of equation (28a) with the boundary condition (30) is 
extremely simple. It is a constant, M = -(l + v)M 0 . 

Nov equation (28b) may be attacked. Introducing the result Just 
obtained, 


S^v 8^w _ Mo 
8x^ 8y^ ^ 


and the boundary condition is, of course, w = 0. 


This differential equation with this boundary condition is known 
in the theory of torsion of a bar having rectangular cross section. All 
that is necessary is to translate the solution known there into the 
terminology of the plate problem. 


The solution is: 


w = 


M o ( y£ 

k'8 2 . 


n-1 cosh 


mtx mty ' 
__ cos - rr *- 


(-D 


1 , 3 > 5 > • • 


3 v nna 
n J cosh 


2b 


(31) 



MCA TN 2612 


57 


One may easily verify that this expression satisfies the differential 
equation and that w = 0 for y = ±b/2. At the other two edges, 
x = ±a/2, equation ( 31 ) yields: 


V = 


% 

K 


1)2 y2 4b 2 

T + T + n 3 





cos 



The sum in the parentheses happens to he the Fourier series representa- 
tion of the function 


fl 2 /b 2 _ y2\ 

4b 2 \ 8 2 / 


valid in the interval -b/2 ^ y ^ b/2, and, therefore, the expression 
in parentheses vanishes at every point of the two edges x = Constant 
of the plate. This proves that the solution really satisfies the 
condition w = 0 on all four sides of the rectangle. 

The bending moments can now be found easily by introducing the 
solution (31) into the elastic law (29) r 


M x = ~Mq 


4(1 - v) ^ 

v + 2 


1,3,... 


(- 1 ) 


n-l , n*x nity 

ii— i cosh -T-— cos — r— 

P b b 


n cosh 


nna 

2b 


My = -M q 


4(1 




These are the bending moments produced in step three. To obtain the 
bending moments for the original problem, the moment Mq of the second 
step must be added; the first step makes no contribution to the stresses. 
When Mq is expressed by the temperature difference T, writing 
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Mq = K(l + V)f 


finally 


M v = 


~ - Eat% 
3« 


Eat% 
12(1 - V) 


00 

f- Z 

1 § 3 9 a • 


2=i cosh 

(- 1) 2 — — 


CQS Hll 


My - 


3« 


Eat^T 


w 

z 

1,3, .. . 


(- 1 ) 


n-1 

2 


cosh cos HI 

D D 


n cosh 


mta 

2b 


The solution would not be complete without having the twisting moment M^y 

Step two does not make- a contribution to it, but it can be obtained from 
equation (3l) alone, using the formula - , 


There is obtained 


V = -K(1 - V) 


d^w 
8x 8y 



— Eat 2 T 
3n 



The formulas for the moments are the solution of the problem as it 
was formulated at the beginning. It is now necessary to discuss this 
result and to draw some practical conclusions from it. 
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At the edges y = +b/2 


n *y A 

cos = 0 

b 


for all odd integers n and hence 


M = — Eat 2 ! 
he 12 


^ = 0 

That the moment My vanishes corresponds to the assumed simple support. 
That M x does not do so, is ^.ue to the fact that the edge is kept straight 
in spite of the applied temperature difference. 

At the edges x = ±a/2 there is the corresponding result 


=0 


— Eat 2n 
12 


My = — Eat^T 


but to obtain it the Fourier series 


00 n-1 

Z / .1 2 l nn; 

(_D _ cos ~= 

1 » 3 > • • • 


must be added up which yields rt/4 for all points of the 
interval -b/2 ^ y ^ b/2. 

' . ■ 

■ 

For interior points of the plate, the series appearing in the 
formulas M x and M y has very good convergence. For the center of 

a square plate, 

M x = My 

= Eat 2 T 
24 


can be found easily. 
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The most interesting part of the solution is the formula for the 
twisting moment M xy . When x < a/2, the quotient of the two hyperbolic 

functions decreases exponentially with increasing n and produces a 
good convergence of the series; but on the edge x = a/2 this beneficial 
influence is lost and the series converges slowly. If x = a/2 is 
kept and the corner y = b/2 is approached, the series becomes 



(- 1 ) 


n-1 
2 1 
n 


n« 
sin -7j- 



and this series is divergent or, if this expression is admitted, yields 
the value » . This singularity oi the twisting moment which, of course, 
appears at all four corners, is of practical importance. It is true 
that real objects always find a way to avoid infinite stresses. Here 
the finite thickness of the plate, the finite width of the zone to which 
the reactions are applied, and the elastic yielding of the support may 
act in this way, but, nevertheless, the singularity in our solution 
reveals the fact that stresses near the corner will be extremely high 
and that it would be wiser to round the corners liberally than to trust 
the good will of the structure. 


Buckling of Cylindrical Panel 

The metal skin of the pressurized cabin is subdivided into rec- 
tangular panels by rings and axial stiffeners (stringers). In every such 
panel, the wall is subjected to a hoop stress a 0 due to the cabin 

pressure p. Additionally, there may be an axial stress cr x (tension 

or compression) and a shear stress t (fig. 4l). The hoop stress a 0 

increases the shear stress t required for buckling in the presence of 
a given cr x . 

For this buckling problem, Kromm (reference 7) has worked out two 
diagrams which give the critical shear t as a function of the hoop 
stress ct 0 assuming either a x = 0 (case "a", fig. 43(a)) or 

°x = a (ZS/^ (case "b", fig. 43(b)). Of course, the ratio between the two 
stresses cj x and may have any other value between or beyond these 

limits, but, since the influence of cr x on the critical shear is not 
large, the choice made by Kromm is sufficient. 

Kromm' s paper gives only a short description of the method used for 
solving the problem, referring for more details to his earlier papers on 



NACA TN 2612 


6l 


stability problems in cylinders. A still shorter outline will be given 
here of the laborious procedure and an explanation of the diagrams 
resulting from it. 

The object of the graphs is a rectangular panel, cut from a circular 
cylinder of radius a and supported on its edges. Its length is supposed 
to be much greater than its width so that the buckling is not influenced 
much by the support on the curved sides. 

When this panel is subjected to the cabin pressure p it will bulge 
out, again forming a cylindrical surface, but with a smaller radius r < a. 
This radius r may easily be found if the lengthwise edges of the panel 
are fixed. In the deformed state (fig. k2) a simple consideration of 
equilibrium yields for the hoop stress the relation 


Od = 




pr 


(32) 


On the other hand, the length of the arc of radius r and chord b is 


l 





the length of the same arc before deformation is 



and hence the hoop strain is 



G i = Ee 0 

- Eh£/L _l\ 

2k \r2 " a 2 J 


When a x = 0 (case "a"). 
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Equating this to the value from equation (32), 

P = tb^ l/l_ _l\ 

E 2k r\r2~ a 2j 

In case "b," E has simply to be replaced by 2E/(2 - v). 

In connection with the buckling problem, the relations between p 
and' r must be expressed by a certain set of dimensionless variables 
used there. Using Kromm' s notation, the curvature of the undeformed 
cylinder is described by the parameter 


a v, 


12(1 - v 2 ) b^ 


A 2 


(33) 


that of the deformed shell, by 


ob = 


12(1 - v 2 ) b 4 


* 4 r-t 


2^2 


and the pressure p, by the parameter 


*p - It - v 2 ) 3/2 4 


(3k) 


The relation just found between p and r reads then in case "a" 


and in case "b": 


^ (<u - «b) = 576 ~ kp 


= 1.035k, 


- (1- 1)576 ^kp 


= 0.880kp 
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After this preparation, the buckling problem for a cylinder of 
radius r may be solved. This may be done by either of the two standard 
methods. The differential equations for the components u, v, and w 
of the displacement may be formulated and solved, or the expression for 
the variation of the potential energy may be established and equated to 
zero for every possible variation u, v, or w of the displacements 
compatible with the boundary conditions. In both cases the displace- 
ments are introduced as double Fourier series, and the buckling condi- 
tion finally assumes the form that a determinant of infinite order must 
be equal to zero. For the numerical evaluation, this determinant is 
approached by a section of moderate size, not necessarily situated at 
its upper left corner. The diagrams in figures 43(a) and 43(b) have 
been computed in this way. 

For the application of these diagrams, it is necessary to know the 
boundary conditions assumed. Since the plate was supposed to be long in 
the x-direction, no conditions were fixed for the curved edges 
x = Constant. On the straight edges, four conditions must be given. 

The following choice was made: Displacement parallel to the edge u = 0, 
radial displacement w = 0, clamping moment (see fig. 8) M0 = 0, and 

additional hoop stress cr^ = 0. 

The first three of these conditions appear to be reasonable at first 
sight. Also the last one is quite usual in buckling problems of this 
kind, but it seems to contradict the assumption of 'unyielding supports 
made for the determination of r. This contradiction may be easily 
resolved. The underlying idea is that the panel is part of the wall of 
a cylindrical fuselage and has many neighbors which are in the same 
situation. When the pressure p is applied to these panels, they will 
all develop the same hoop stress cr^. Although the stringers usually 
have but little bending stiffness, they cannot deflect in the v-direction 
(fig. 4l), because the forces applied to them from both sides are in 

equilibrium. When buckling occurs, the situation is quite different. A 
system of folds is formed' in each panel, and with them additional hoop 
stresses are set up. If the stringer were very stiff, corresponding 

forces ajt would be transmitted to it, and they would not be the same 
from both sides but would pull at some places to one side and at other 
places to the other side. When the stringer is weak, as is usual, it 
will deflect so much that the stresses become almost zero, and the 
safest assumption for the determination of the buckling stress is the 
one made . 

The use of the diagrams may now be described. From the given data 

and kp are computed according to equations (33) and (34) and 

located in the diagram. The values kp are found at the left side of 
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the graph and refer to the more-or-less horizontal curves. For c%, 
some values are given along the upper edge. The curves to which they 
refer are almost vertical and join at their lover end, more or less 
tangentially, a vertical line bearing the same number on the \f u> scale. 
From this relation it is easy to interpolate more curves ^ and to 
find the values for those curves which are not numbered in the 

graph. It also appears that in the right-hand half of the diagram the 
\fa ^ curves are practically identical with the vertical coordinate lines. 


From the point which corresponds to the given values of kp and 
cu 0 follow a horizontal line toward the left and read there t/ct*. 


When multiplied by the reference stress 


a* = * 


E 


3 1 _ V 2 b 2 


it yields the critical value of the shear stress t. 


If, incidentally, cr x corresponds to one of the two cases "a" 

or "b," it is necessary to consult only one of the diagrams. For 
other values of a x it is necessary to use both and then find the final 
value t by interpolation. Since the influence of cr x is not great, 
even an extrapolation will be possible, within moderate limits. 


Stanford University 

Stanford, Calif., September 7, 1950 
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Figure 3.- Part of cross section through a cylindrical shell with stringers 


( 



Figure iu- Factor k for longitudinal stress in skin of stiffened shell. 
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(b) Rings far apart, light stringers. 

Figure 10.- Typical distribution of bending moment and hoop force H0 

in shell between two rings. 










NACA TN 2612 




N«p 



(b) Cylinder, bulkhead, and ring taken apart. 

Figure 17 .- Axial section through a cylindrical shell and a bulkhead. 
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'Figure 18. -^Meridional force N 0 and hoop force Nq in two different 
bulkheads. Upper half, extremely flat shape 5 lower half, sufficiently 
dished shape. 
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(a) Parts before pressure is applied* 



(b) Exaggerated scheme of deformation when pressure p is applied to 

cabin at left of bulkhead* 

Figure 21.- Continuous cylindrical shell having _ an inserted bulkhead. 
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Figure 29.- Sign convention for forces transmitted from spherical shells 

to connecting ring. 



Figure 30 .- Element of ring shown in figure 29. 
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Figure. 31.- Distribution of circumferential force Nq and normal force 
at. edge of ellipsoidal bulkhead shown in left half of figure 26. 











leflection of 
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to stiffener 










Figure UO.- Rectangular plate and plate element 







12C 


NACA TN 2612 


89 


\ ^ \ v.\ \ 



m m 


Figure Ul.- Cylindrical panel. 



Figure I 4 . 2 . — Section through cylindrical panel before and after application 

of internal pressure p„ 
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in figure 41. 
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NACA - Langley Field, Va. 


Figure 1*3.- Concluded 



